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We study the expansion dynamics of a quark-antiquark plasma droplet 
from an initial state with restored chiral symmetry. The calculations are 
made within the linear a model scaled with an additional scalar field repre- 
senting the gluon condensate. We solve numerically the classical equations 
of motion for the meson fields coupled to the fluid-dynamical equations for 
the plasma. Strong space-time oscillations of the meson fields are observed 
in the course of the chiral transition. A new phenomenon, the formation and 
collapse of vacuum bubbles, is predicted. The particle production due to the 
bremsstrahlung of the meson fields is estimated. 

Introduction: It is commonly believed that the conditions for chiral symmetry restoration 
and colour deconfinement can be reached in the course of an ultra-relativistic heavy-ion 
collision. The Quark-Gluon Plasma (QGP) is expected to be formed at some intermediate 
stage of the reaction. Since strong collective expansion may develop already in the QGP, its 
subsequent transition to the hadronic phase should be treated dynamically [|TJ. At present 
this is possible only on the basis of effective models obeying the symmetry properties of 
QCD. 

Our considerations below are based on the linear a model which respects approximate 
chiral symmetry. In addition to the usual chiral fields, a and iz = {ni, tt2, ^3}, the model 
includes the dilaton or glueball field, \i to simulate the trace anomaly of QCD 0. The a 
and x represent the quark and gluon condensates as effective meson fields. Models of this 
kind were used earlier for nuclear matter (see e.g. |§). 



The dynamical model: The effective chiral Lagrangian for constituent quarks interacting 
with the background meson fields is written as 
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Here U is the usual Mexican Hat potential scaled by the glueball field x (below we take 
n = 3) and W is the effective potential responsible for the scale invariance breaking. The 
trace of the energy-momentum tensor for the above Lagrangian, Tjf = — A 2 x 4 , is assumed 
to be proportional to the gluon condensate, {G 2 ). The scale parameter A is of the order of 
A Q cd ~ 200 MeV. 

This Lagrangian leads to the normal vacuum state where chiral symmetry is spon- 
taneously broken: cr = f n = 93 MeV, 7r=0, X = Xo = 136 MeV. The parameters 
of the Lagrangian are chosen so that in the normal vacuum the constituent quark mass 
m q — 9 fir — 313MeV, the a-meson mass m 2 = 2A 2 / 2 + m 2 « (0.6 GeV) 2 and the glueball 
mass m 2 G = ^X^Xo + 0{ m V) ~ (1-7 GeV) 2 . The energy density associated with breaking the 
gluon condensate, B = X\xV^ is fixed to 0.5 GeV/fm 3 Q. In the case of thermal equilib- 
rium and a frozen x field, x = Xo-, the model leads to a chiral transition at temperatures 
around 130 MeV §. 

The form of the effective glueball potential eq. flU) is motivated also by the instanton 
liquid model (see the recent review ||) if x 4 is identified with the instanton density. As pre- 
dicted by this model, the instanton density is significantly suppressed at high temperatures. 
In our calculations the coupling of the gluon condensate to the thermal bath is parametrized 



in a simple form J(x)x 2 {x) (the last term in eq. @), where J{x) = AT 2 (x) and T(x) is the 
local temperature. The coupling strength A w 2.4 is chosen so that the gluon condensate, 
X 4 , is reduced by about 60% at T = 280 MeV. 

The model presented above is not fully consistent since it does not include exsplicitely 
the thermal gluon excitations. In doing so we were motivated by the recent analysis || of 
lattice data demonstrating that at temperatures of interest here the gluons have a rather 
large effective mass of about 0.6-0.8 GeV. Therefore, the contribution of thermal gluons to 
all thermodynamical quantities is significantly suppressed by the Boltzman factor. On the 
other hand, their contribution is included implicitely in the coupling term Jx 2 determining 
the degree of reduction of the gluon condensate at high temperatures. 

Below we adopt the mean field approximation considering a, ir and x as classical fields. 
The equations of motion for these fields are obtained by applying the variational principle 
to the above Lagrangian. The source terms in these equations are determined by the dis- 
tribution of quarks and antiquarks, which in principle should be found by solving the Dirac 
equation. Due to the interaction with meson fields, quarks acquire an effective mass J7| 

m q (x) = g\Ja 2 (x) + tt 2 (x), (3) 

which, in general, is space and time dependent. This makes an exact solution of the Dirac 
equation very difficult. To avoid this problem one should make further approximations. 
A reasonable starting point is the Vlasov-type kinetic equation for the scalar part of the 
quark- ant iquark Wigner function f(x,p), 
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where I co u is the collision term. In refs. Q,§,|TIj] the collisionless (I co ii — 0) approximation was 



used for the propagation of quarks in background meson fields. Obviously this approximation 
can be justified only for the late stages of the expansion. 

Here we consider another approximation which is more appropriate to high temperatures. 
Namely, we assume that the partonic collisions are frequent enough to maintain local ther- 



modynamical equilibrium. In this case f(x,p) can be represented in terms of the equilibrium 
distribution functions characterized by local temperature T(x) and chemical potential 

By multiplying both sides of eq. (^) with p^, projecting on the mass shell , p^p^ = m 2 (x), 
and integrating over 4-momenta one arrives at the equations of relativistic hydrodynamics 
(see also JT2]]) 
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where T^ v [x) is the energy-momentum tensor, T^ v = (S + P)w M M iy — Pg^ u , and u v (x) is 
the collective 4- velocity of the quark- ant iquark fluid. Here the energy density S, pressure 
P and scalar density p s are functions of T(x), and m q (x). They can be expressed in 
the standard way through the Fermi-Dirac occupation numbers. We solve these equations 
consistently with the equations of motion for the meson fields which determine the quark 
effective mass through eq. (^|). This is why we call this approach Chiral Fluid Dynamics 
(CFD). 

The evolution of the glueball field is driven by the couplings to the chiral fields and to 
the thermal bath. The corresponding source J{x) drops with the characteristic hydrody- 
namical time of order of a few fm/c. Due to some uncertainties in the glueball Lagrangian, 
particularly, in the derivative terms, below we consider two options of the model , i.e. full 
dynamics as described above and pure chiral dynamics with a frozen glueball field (x = Xo)- 

Numerical results: For numerical simulations we have used the RHLLE algorithm de- 



scribed and tested in JT3[ for fluid dynamics and the staggered leapfrog method for the 
field equations. Below we present results for the real-time evolution of spherical droplets 
of radii R = 2 fm and R = 4 fm. In the initial state we take a baryon-free (/i = 0) fluid 
with a Woods-Saxon temperature profile and a linear profile of the collective momentum 
density. Initially the system is assumed to be in the chiral- symmetric phase at temperature 
T « 280 MeV. The initial conditions for the fields are chosen uniformly within the droplet 
and smoothly interpolated to their vacuum values outside the droplet. We assume that a 
and x fields are initially close to their equilibrium values at this high temperature. 



The results of the calculations for R = 4 fm are presented in Figs. 1 and 2. It is seen 
that within a few fm/c the energy density of the fluid drops from the initial value of about 
5.0 GeV/fm 3 to below 0.1 GeV/fm 3 . A shell-like structure of the matter distribution is 
clearly seen at late times P^JI5| . The fluid is cooled down to T = 130 MeV already at 
t ~5 fm/c. As Fig. 2 shows, at this time the a field changes rapidly from its initial value, 
a ~ 0, towards the new asymptotic value, o = f n . This transition is accompanied by strong 
nonlinear oscillations of the all coupled fields. 

The pion field oscillations are especially strong and spread over the whole space within 
the light cone. In accordance with previous studies [|16H 20[], our calculations show that soft 
pion modes are indeed strongly amplified (by a factor 10) in the course of the chiral transition 
even in a finite expanding droplet. As suggested earlier (see e.g. |21~|-f2"3|j ) , a perfect isospin 



alignment of the classical pion field should lead to a non-statistical distribution of the ratio 
of neutral to charged pions. One should however bear in mind that these coherent pions 
will be accompanied by a large number of genuine pions (in the considered example, about 
1000) produced at the hadronization of plasma. 

It is interesting to note that the heavy a and x fields have quite different dynamics com- 
pared to the pion field. Initially they evolve almost adiabatically following the instantaneous 
temperature. Instead of expanding they first shrink and then rebound at about the time 
of the chiral transition, when strong oscillations start. One can understand this behaviour 
from the following consideration. The regions where a and x significantly deviate from their 
vacuum values, f n and Xo? have an extra energy density A£ vac ~ B, compared to the normal 
vacuum. This vacuum energy excess generates a negative pressure, P vac = —A£ vac . Such 
regions can survive only untill the internal pressure of matter (in our case, quark- ant iquark 
fluid), P m ati is large enough to counterbalance the external vacuum pressure, i.e. when 
Pmat + Pvac > 0. This condition is always fulfilled in an equilibrated system, Howevere, 
in the course of a rapid expansion the opposite condition, P mat + P vac < 0, can eventually 
be reached in a certain region of space which we call a vacuum bubble. Then the outside 
vacuum propagates into this bubble trying to minimize its size. This process looks like a 



collapse of an air bubble in a liquid. In our case the role of a liquid is played by the vacuum 
quark and gluon condensates. 

The collapse starts from the surface of the quark-antiquark droplet. As the energy 
density of the fluid decreases the speed of the ingoing wave increases. Finally the true 
vacuum penetrates to the center. Due to the inertial forces the condensates overshoot their 
equilibrium vacuum values. This is why very strong oscillations are developed at the center 
of the bubble. This violent dynamics may lead to very interesting phenomena like particle 
production by the bremsstrahlung mechanism, reheating of the fluid or trapping of some 
quarks and antiquarks (especially heavy ones s, s, c, c) in the bubble. In this paper we 
consider only the first process. 

Particle production: In general, a time-dependent meson field 0(r, t) can be represented 
asymptotically as an ensemble of quanta of this field. Using the coherent state formalism 
one can write the explicit expression for the momentum distribution of the mesons produced 
(see e.g. refs. [0,|TJ) 
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where = J k 2 + is the single-particle energy of a meson with the vacuum mass m^. In 
this formula 0(k, t) and 0(k, t) are the 3-dimensional Fourier transforms of the meson field 
4>(r,t) and its time derivative <j)(r,t). They are obtained from the dynamical simulations 
described above. The r.h.s. should be calculated at sufficiently late times when nonlinearities 
in the field equations are negligible. 

The numbers of produced particles depends sensitively on the initial conditions for the 
fields, droplet size and the expansion dynamics. In Table 1 we present the particle numbers 
(calculated at t = 30 fm/c) for the illustrative example of Fig. 2, R = 4 fm, with the initial 
glueball field 0.79xo (case I) and O.lxo (case II), as well as for a smaller droplet, R = 2 
fm. For comparison, the results for the frozen glueball field, \ = xo, are also shown. The 
inclusion of a dynamical glueball field leads to an increase in the pion number by a factor of 
2 or more, mainly due to the additional amplification of the "pion laser" modes fH|. They 

6 



originate from the splitting of soft a modes into pion modes with momenta of about m a /2. 
But a larger fraction of energy, associated with the initially suppressed gluon condensate, 
goes into the bremsstrahlung of a and x fields. One can see that much more particles, 
especially glueballs, are produced when the system is initially trapped in a metastable state 
with x ~ (case II). On the other hand, significantly fewer particles are produced from a 
smaller droplet with R = 2 fm. 

Conclusions: It is demonstrated that the chiral transition in an expanding finite droplet 
is accompanied by strong space-time oscillations of the background fields. The long wave- 
length modes of the pion field are strongly amplified, by a factor 10, in the course of transi- 
tion. The gluon condensate brings into play a new scale of energy density, B = 0.5GeV/fm 3 , 
which changes significantly the dynamics. 

The simulations reveal a novel phenomenon, the formation and collapse of the vacuum 
bubbles, associated with the regions of out-of-equilibrium quark and gluon condensates 
and low matter pressure. The additional energy released in the collapse goes partly into 
the coherent pions, but to a larger extent, to the production of a- mesons and glueballs. 
Due to a very large width of the a meson its direct observation in heavy ion collisions 
is practically impossible. But the glueballs can be detected by the characteristic decay 
channels G — > nn, KK with widths of about 100 MeV, as well as by the electromagnetic 
decay (G — > 77) with a width of a few keV. 

In the future we are planning to improve the model in two directions. First, a more 
realistic treatment of the gluon condensate and its coupling to the partonic plasma should 
be introduced. Second, friction terms due to the interaction of meson fields with the thermal 
bath should be included. 
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TABLE I. The numbers of neutral pions {N no ), a mesons (N a ) and glueballs (Nq) produced 
by the brensstrahlung mechanism in the course of a quark-antiquark droplet expansion. Results 
are presented for initial droplet radii 2 and 4 fm, as well as for full (cases I and II) and frozen \ 
dynamics. 
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FIG. 1. Contour plots in the r — t plane (case I) of the energy density in GeV/fm 3 (top from 
left the subsequent decreasing levels are: 5.0,4.55,4.05,3.55,3.05,1.55,2.05,2.55,1.05,0.55,0.35,0.15 
and 0.05 GeV/fm 3 ) and temperature in MeV (bottom, from left the subsequent decreasing lev- 
els are: 280,260,240,220,200,180,160,140,120,100,80,60 and 40 MeV) in the expanding spherical 
quark-antiquark droplet of initial radius R = 4 fm. The initial collective velocity is about 0.2c at 
the surface. 
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FIG. 2. The evolution of pion, sigma and glueball fields in the t — r plane for an expanding 
droplet of initial radius R = 4 fm. The initial values ((case I) of the fields are: a = —0.05f w , & = 0, 
7T 3 = 0.2f n , ir 3 = 0,x = 0.79xo and * = 0. 
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